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ABSTRACT 

An experimental test facility Is used to measure the leakage and 
rotordynamic coefficients of constant**c 1 earance and convergent- 
tapered annular gas seals. The results are presented along wltH 

i 

the theoretically predicted values. Of particular Interest is the 
prediction that optimally tapered seals will have significantly 
larger direct stiffness than straight seals. The experimental 
results verify this prediction. Generally the theory does quite 
well, but falls to predict the large Increase In direct stiffness 
when the fluid Is prerotated. 
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NOMENCLATURE 


A 

C r 

CtC 

F 

f 

K, k 

k l 

M 
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mo, no 


P 

R 

R 


u z ,u e 

V m 


Y 


y 
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T 

to 


cross-sectional flow area of the seal (L 2 ) 

seal clearance (L) 

■ 

direct and cross-coupled damping coefficients (FT/L) 
reaction force (F) 

Darcy friction factor 

direct and cross-coupled stiffness coefficients (F/L) 

loss coefficient 

Mach number 

mass flow rate (M/T) 

Hlrs' constants 
2 

pressure (F/L ) 
shaft radius (L) 

Reynolds number, 2pV m c p /y 
perfect gas constant (L 2 /T 2 ©) 

fluid velocity In the Z and 0 direction (L/T) 

bulk-flow velocity of the fluid (L/T) 

specific heat rato 

viscosity (FT/L 2 ) 

density (M/L 3 ) 

shear stress (F/L 2 ) 

shaft angular velocity (T - *) 


Subscripts 

a,b reservoir and sump 

s,t = straight and tapered 

X,Y,Z = rectangular coordinate directions 


INTRODUCTION 


Part 1 of this paper described the deslon, development, and 
operation of a test apparatus and facility which have been 
developed to measure leakage and rotordynamlc coefficients of 

annular gas seals. This test apparatus has been designed and used 
to measure, rotordynamlc coefficients of plain annular seals (as 

used In f ioatlng-rlng seals), plain seals with honeycomb stators, 
and labyrinth seals. This part (Part 2) of the paper provides a 
comparison of test results and theory for plain annular seals 
having constant-clearance and convergent-tapered geometries. 
Subsequent publications will provide test results for honeycomb 
and labyrinth seals. 

As described In Part 1, the rotordynamlc coefficients for a 
gas seal are defined by the following linearized force- 
displacement model. 



Where (X, Y) define the relative motion of the seal rotor relative 

to Its stator, (F , F ) are the components of the reaction force 

X Y 

acting on the rotor, and (K f K , K , K ) and (C , C , C , 

XX YY XY YX XX YY XY 

C ) are the stiffness and damping coefficients respectively. 
YX 

Equation (1) applies for small motion of the rotor about an 
arbitrary eccentricity position. For small motion about a 

centered position, the following simpler model applies. 



Although the tost facility has the capability of measuring the 
eccentric-position rotordynamlc coefficients of Eq. (1), present 
analytical solutions are only capable of predicting the 
centered-pos It Ion coefficients of Eq. (2). Thus, all test results 
reported In this paper will be for the shaft-centered position. 
Analytical Developments 

Literature pertaining to the analytical solution for the 
rotordynamlc coefficients of annular gas-path seals Is extremely 
limited. Fleming [1,2] analyzed convergent-tapered gas seals 
using a seperate analysis to determine the direct stiffness K and 
the direct damping C. His analysis concludes that optimally 
tapered seals will develop considerably higher direct stiffness 
than straight seals. However, since his analysis does not Include 
fluid prerotation (preswirl) he was not able to predict the 
cross-coupled coefficients c and k. Recently Nelson [3,4] 
developed an ana 1 yt I ca 1 -computat 1 ona 1 method to solve for all of 
the rotordynamlc coefficients In Eq. 2. solution technique 

Is similar to that developed by Childs [5»6]. That Is, a set 
of governing equations based on Hlr's [7] turbulent bulk flow 
model are developed, and then a perturbation analysis Is employed 
to obtain a set of zero- and first-order equations. Integration 
of the zeroth-order equations yields the leakage, and Integration 
of the first-order equations yields the direct and cross-coupled 
rotordynamlc coefficients. The analysis handles choked or 
unchoked flow, the effect of different stator and rotor 



surface roughness conditions* and straight or tapered seal 
geometries. His solution agrees with Fleming's prediction that 
convergent -tapered seals develop higher direct stiffnesses. The 
test results presented In this paper are compared with the 
theoretical predictions of Nelson's solution. 



TEST VARIABLES 


There are fifteen Independent Input variables required for 
the theoretical solution of the smooth* convergent-tapered, 
annular gas seal. These values can be grouped Into the three 
categories: geometry, gas properties, and operating conditions. 

The values (or range of values) used for the theoretical solution 

are shown in Table I below. Following the table is a description 
of how these values were obtained from the test apparatus. 

Geometry 

length: 8 cm (2.00 In) 

diameter: 24 cm (6.00 In) 

entrance clearance: 0.7366 mm (29 mil) 


exit clearance: 0.7366 mm (29 mil) 

1 . 1 1 43 mm ( 45 mi 1 ) 


( stra 1 ght ) 
(tapered) 


Hlr's constants: no = 0.187 mo *= -0.333 

Gas Properties 

viscosity: 1.9xl0 -5 N*s/m 2 (4.0x10-7 lb-s/ft Z ) 

gas constant: 287 m»N/kg»K (53.3 ft* lb/lb m *°R) 

specific heat ratio: 1.4 

Operating Conditions 

reservoir pressure: 1.7 to 7.2 bar (25 - 104 psl) 
reservoir temp: 305 K (90 °F) 


sump pressure: 
preswl r 1 ratio, U eo (0)/Rw: 
shaft speed: 
process Iona 1 speed* 


1 .0 bar (14.7 psl) 
-81 to +81 
200 to 8000 cpm 
0.98 and 1 .24 cpm 


Table 1. Values of the Independent Input Variables 



Geometry, Figure l shows a simplified cross-sectional view of the 
rotor and stator. Both the stator and rotor are fabricated from 
304 stainless steel with a 4 win finish. The seal's length* 
diameter* and clearances are obtained from the given dimensions. 

The theoretical solution Is oased on Hlrs' turbulent bulk- 
flow model. That Is, In the relationship for the wall shear stress, 

x * <l/2)pV m 2 f (3) 

.It Is assumed that the Fanning friction factor, f, Is a variable 
according to the following relationship: 

f * no(R) mo (4 ) 

In these equations^ V m is the mean bulk-flow velocity and R is the 
Reynolds number, both relative to the surface upon which the shear 
stress Is acting. The Hlrs' coefficients no and mo account for 
the surface roughness and must be empirically determined. To 
determine the Hlrs' constants for the test aparatus, static tests 
(l.e. no shaking nor rotation of the rotor) were run with the 
constant-clearance seal at several different reservoir pressure 
levels. The data collected from these tests Included the mass 
flow rate, m, reservoir temperature, Ta, reservoir pressure, pa, 
and the pressures at five axial locations. From these values, the 
Reynolds number Is calculated, and the Mach number, M, Is 
determined at each of the five axial locations from the following 
equat I on . 


m t= pAM Jy[ 1 + 
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These mach numbers were then used to solve for the left-hand side 
of the standard equation for choked » adiabatic, compressible flow 
In a duct with friction. 

2fL* = + (2+1) W _„Lytl)H 2 _ \ ( 6 ) 

(T yM^ 2y }2[1 + (v-1 )h z 1 J 
• 2 
Here, L Is the length of seal from the axial location m to the 

exit which would be necessary to produce choked flow. The average 

friction factor between any two arbitrary locations m and n can 

then be found from 

« « 

f = C 

where L _ Is the distance between points m and n. The mean 

ill f 1 1 

friction factor was then found by averaging the results of 
different point combinations. Finally, the Hirs' constants, no 
and mo, were evaluated from a power law curve fit of the mean 
friction factors and the Reynolds number obtained at the various 
reservoir pressure levels. 


Gas Properties? The fluid passing through the seal was air. The 
gas properties used were taken from standard tables. 


Operat 1 np 

Conditions: 

To 

determine the 

effect 

of 

fluid 

prerotation 

(preswi r 1 ) , 

the 

test rig design 

al 1 ows 

for 

the 


Installation of a set of Inlet vanes which accelerates the flow In 
the circumferential and axial direction. For tests In which there 
was no prerotation, no guide vanes were Installed. With no vanes, 
the Inlet cross-sectional area Is significantly larger that the 
sealed area, and thus the measured Inlet pressure and temperature 
fs essentially equal to stagnation (reservoir) pressure and 
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( 7 ) 



temperature* and these values were used directly In the 
theoretical solution. For the case In which the guide vanes were 
Installed* appropriate corrections were made to the accelerated 
flow to convert the measured values back to the reservoir 
conditions. 


preswlrl ratio. 

U (0)/Ru, 

1 s 

the 

ratio of 

the entering 

circumferential 

velocity 

to 

the 

surface 

speed of the 


shaft. Although the Inlet guide vane angle Is fixed, this ijatlo 
varies as the pressure drop across the seal and the shaft speed Is 
changed. For the largest pressure drop and the slowest shaft 
speed, the ratio has a value of a±81. The plus sign Indicates 
flow which has been preroteted In the direction of shaft rotation, 
end the negative sign Indicates prerotation against. For the 
smallest pressure drop and the fasted shaft speed, the ratio has a 
value of “±1.8. 

The theoretical model assumes that the shaft precesses In a 
small elliptical orbit about the shaft-centered position. In the 
test apparatus, the shaft Is shaken horizontally about the 
centered position. The shake frequencies used were 58.8 and 74.4 
Hertz (0.98 and 1.24 cpm) . Text results were essentially the same 
for both these frequencies. 



TEST RESULTS 


Figure 2 shows the axial pressure gradient for the static 
test of the constant-cl ance seal. The upper plot Is for 
non-preroteted fluid, and the lower plot Is for prerotated fluid 
In which the prerotation opposes the rotor rotation. The solid 
line Is for the measured values, and the dashed line represents 

theory. For non-prerotation, agreement between theory and 
experiment Is excellent. This Is to be expected, however, since 
the Hlrs' coefficients used In the analysis come directly from the 
measured pressures. For the prerotated case, the agreement Is not 
as good. As stated earlier, when the Inlet guide vanes are 
Installed, appropriate corrections are made to the accelerated 
flow to convert the measured inlet flow conditions back to their 
reservoir conditions. Since the experimental plot shows the 
measured static pressure (not yet corrected to reservoir 
conditions), there Is naturally a difference between theory and 
experiment at the front of the seal entrance. Just after the seal 
entrance, however, theory and experiment should agree. The sudden 
drop In the pressure at the entrance Is due to conversion of 
pressure energy to kinetic energy and due to non- 1 sent ropl c 
losses. The non-1 sentropl c losses are accounted for In the 
theoretical solution by a loss coefficient which Is approximated 
from the data of Delssler [8] as 


1 


1.0 


k , = J 5. 3/ log R 

1 I 10 


( 8 ) 



From the lower plot In Fig. 2, It appears that this coefficient 
underestimates the entrance losses when the fluid Is prerotated 
with the Inlet guide vanes. This same general result Is seen for 
prerotetlon In either direction# and for both the constant- 
clearance and convergent-tapered seal. 

Experimental data was taken at shaft rotor speeds of 200, 

500, and 1000-0000 rpm. In 1000 rpm Increments. Figure 3 shows 
the direct stiffness, Kxx, vs. these rotor speeds for a fixed 
pressure ratio (pa/pb). It can be seen that there Is very little 
dependence on rotor speed. This result Is typical for all the 
rotordynamlc coeff Ichwcs and leakage at any of the pressure 
ratios tested and for both the constant-clearance and convergent 
tapered seal. 

Figures (4-8) shows the measured and predicted leakage and 

rotordynamlc coefficients vs. pressure ratio. In each figure, the 

upper plot is for the constant-clearance seal and the lower plot 

Is for the convergent-tapered seal. Each plot shows the results 

for three cases: (1) no fluid pre>‘otat Ion, (2) fluid prerotated 

In the direction of the shaft rotation, and (3) fluid prerotated 

against the direction of the shaft rotation. Since there was 

little dependence of the leakage and rotordynamlc coefficients on 

shaft speed over the range tested, average values measured over 

the speed range were used as the data points on the plots, with a 

# 

vertical line at each point Indicating the variance. The predicted 
values are shown by the solid line. 

Predicted and measured values for the leakage of the straight 
and tapered seal are shown In Fig. 4. The agreement Is excellent. 
Figure 5 shows the direct stiffness. The most notable result In 
this figure Is the effect that prerotation of the fluid has on the 



experimentally measured values of K. Prerotatlng the fluid causes 
a two- to three-fold Increase* However, the theoretical solution 
does not account for this effect. The predicted values are all 
within the range of the measured values, but essentially Indicates 
very little effect of prerotation. 

Figure 6 shows the cross-coupled stiffness. Agreement 
between the predicted values and the measured values Is fairly 
good. For no prerotation, both measured and predicted are 
nominally zero. For the prerotated, the predicted values run 0 - 
40% larger than measured. 

Figure 7 shows the direct damping. For the straight seat the 
agreement between theory and experiment Is excellent. For the 
tapered seal the theory runs 0 - 30% smeller than measured. 

Figure 8 show the cross-coupled damping. Again, for no 
prerotation, both measured and predicted values are nominally 
zero. However, the theory significantly underestimates c for most 
of the prerotated cases. 



DISCUSSION 


Ab stated In the Introduct Ion, both Fleming's [2,3] and 
Nelson's [4,5] analytical solution predicts that tin optimally 
tapered seal will develop considerably higher direct stiffness 

than straight seals. The experimental results verifies this 
prediction. This can more easily be seen In Fig, 9. In this plot 
the vertical axis represents the ratio of the direct stiffness 
for the tapered seal to the direct stiffness of the 

constant-clearance seal. The theory predicts stiffness ratios of 
1.3 to 2.0 as the pressure ratio Is Increased, and Is essentially 
not Influenced by fluid prerotation. The measured results show 
ratios of 1.2 to 1.5 and 1.6 to 2.3 for prerotated fluid, and 
of 1.6 to 2.7 for no prerotation. 

Generally the theory appears to work quite well. It does an 
excellent Job of predicting the leakage and the direct damping, 
and does a fairly good Job of predicting the cross-coupled 


st 1 ffness . 

The 

most 

notab 1 e 

fa 1 1 ure 

Is In accounting for 

the 

effect of 

fluid 

prerotation 

on the 

direct stiffness. 

As 

prev l ous 1 y 

noted. 

the 

emp l r 1 ca 1 

mode 1 of 

Delssler's [9] used 

for 


the Inlet losses In the theoretical solution Is not adequate when 
the Inlet guide vanes are Installed. Improved modeling of the 
Inlet losses may Improve the accuracy of the stiffness 
predl ct Ions. 

Finally, agreement between experiment and theory for the 
cross-coupled damping coefficient Is at best marginal. However, 
the actual force produced by the effect of c Is much smaller than 
the other rotordynamlc coefficients. Thus, Its effect on 
rotordynamlc calculations Is rather Insignificant. 
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